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ABSTRACT 

We present first results from the 3D-HST program, a near-IR spectroscopic survey performed with the Wide 
Field Camera 3 on the Hubble Space Telescope. We have used 3D-HST spectra to measure redshifts and Ha 
equivalent widths (EWh q ) for a complete, stellar mass-limited sample of 34 galaxies at 1 < z < 1 .5 with M star > 
10 11 M Q in the COSMOS, GOODS, and AEGIS fields. We find that a substantial fraction of massive galaxies 
at this epoch are forming stars at a high rate: the fraction of galaxies with EWh > 10 A is 59%, compared 
to 10 % among SDSS galaxies of similar masses at z = 0.1. Galaxies with weak Ha emission show absorption 
lines typical of 2-4 Gyr old stellar populations. The structural parameters of the galaxies, derived from the 
associated WFC3 F140W imaging data, correlate with the presence of Ha; quiescent galaxies are compact 
with high Sersic index and high inferred velocity dispersion, whereas star-forming galaxies are typically large 
two-armed spiral galaxies, with low Sersic index. Some of these star forming galaxies might be progenitors 
of the most massive SO and Sa galaxies. Our results challenge the idea that galaxies at fixed mass form a 
homogeneous population with small scatter in their properties. Instead we find that massive galaxies form a 
highly diverse population at z > 1, in marked contrast to the local Universe. 
Subject headings: cosmology: observations — galaxies: evolution 
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1. INTRODUCTION 

In the nearby Universe galaxies with stellar masses > 
10 n M Q form a homogeneous population, with small scat- 
ter in their properties at fixed mass (e.g., Djorgovski & Davis 
1987; Blanton et al. 2003; Kauffmann et al. 2003a). This 
homogeneity is somewhat puzzling in the context of standard 
models of galaxy formation, as in these models star formation 
and merging extend to recent epochs (e.g., De Lucia et al. 
2006). Models that aim to explain the tight scaling relations 
of massive galaxies usually invoke "dry" mergers in combina- 
tion with feedback from AGN to prevent gas cooling and star 
formation (e.g., Croton et al. 2006). 

Many studies have measured the properties of massive 
galaxies at earlier cosmic epochs to better constrain when and 
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how they were formed. Interestingly, it was found that mas- 
sive galaxies have red rest-frame optical colors at least out to 
Z ~ 2 (Bell et al. 2004; Faber et al. 2007; Brammer et al. 
2011) and probably beyond (e.g., van Dokkum et al. 2006; 
Marchesini et al. 2010). Many of these galaxies are com- 
pact quiescent galaxies with little or no ongoing star forma- 
tion (e.g., Daddi et al. 2005; van Dokkum et al. 2008; Dam- 
janov et al. 2009), but at redshifts z > 1 a subset of the popula- 
tion has MIPS 24 /im fluxes and near-IR colors characteristic 
of dust-obscured star-formation (e.g., Papovich et al. 2006; 
Williams et al. 2010; Brammer et al. 2011). These galaxies 
mascarade as "dead" red sequence galaxies as they have very 
similar optical colors but in fact have very high inferred star 
formation rates. 

These apparently star-forming, massive galaxies are pro- 
genitors of at least a subset of massive galaxies today. In 
the nearby Universe strongly star forming galaxies are typ- 
ically gas-rich mergers, but several studies have argued that 
at higher redshift such galaxies more resemble "scaled-up" 
spiral galaxies than local mergers (Wolf et al. 2005; Muzzin 
et al. 2010; Elbaz et al. 2011). Interpreting these galaxies, 
and their relation to the compact quiescent galaxies which ex- 
ist at the same epoch, has been hampered by a lack of mass- 
complete samples with homogeneous data at the redshifts of 
interest. Most existing samples are luminosity- rather than 
mass-selected, are based on photometric redshifts, lack rest- 
frame optical morphological information, and/or lack well- 
calibrated star formation diagnostics. 

In this Letter, we construct and study a spectroscopic stel- 
lar mass-limited sample of galaxies at 1 < z < 1.5, in order 
to quantify the properties of massive galaxies when the Uni- 
verse was k 5 Gyr old. We use data from the 3D-HST survey 
(GO- 12177, GO- 12328), an HST/WFC3 Treasury program 
that provides rest-frame optical spectroscopy and imaging at 
0"13 resolution for thousands of distant galaxies. 
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2. DATA AND ANALYSIS 

The data presented here were obtained in the context of 3D- 
HST, a 248-orbit HST Treasury program in Cycles 18 and 
19. 3D-HST is a wide-field, two-orbit depth survey of four 
(GOODS-South, UDS, AEGIS, and COSMOS) of the CAN- 
DELS fields (Grogin et al. 2011, Koekemoer et al. 2011) 
with the Wide Field Camera 3 (WFC3) G141 grism. The 
fifth CANDELS field, GOODS-North, was observed by the 
Cycle 17 program GO-11600 (PI: Weiner); these data are in- 
cluded in the 3D-HST project. The grism provides spatially- 
resolved spectra of all objects in the WFC3 field, covering 
the wavelength range 1.15 /im - 1.65 /im with point-source 
spectral resolution R w 130. Accompanying direct imaging 
in the F140W filter provides the necessary information for 
wavelength calibration, as well as photometric and morpho- 
logical information. A typical two-orbit pointing comprises 
« 5100 s of integration time in the G141 grism and w 800 s 
in the F140W direct imaging filter. In addition to the WFC3 
near-IR spectroscopy and imaging, 3D-HST provides parallel 
ACS G800L grism spectroscopy and accompanying F814W 
imaging in the optical. The survey is described in G. Bram- 
mer et al., in preparation. 

The F140W data were reduced using standard procedures 
for WFC3 imaging data, and resampled to a rectangular out- 
put grid with 0"06 pixels. This fine sampling is justified as 
the four individual exposures in each visit are shifted with re- 
spect to each other by an integer + 0.5 pixels in x and y. Spec- 
tra were extracted from the grism data using the aXe soft- 
ware (Kiimmel et al. 2009), in a similar way as described 
by van Dokkum & Brammer (2010) and Atek et al. (2010). 
Details on the specific 3D-HST reduction will be provided in 
G. Brammer et al., in preparation. The spectra were combined 
with photometry at other wavelengths, using publicly avail- 
able photometric catalogs of Whitaker et al. (2011) (COS- 
MOS, AEGIS), Kajisawa et al. (2011) (GOODS-North), and 
Wuyts et al. (2008) (GOODS-South). Next, redshifts and 
emission line fluxes were determined from the combined pho- 
tometric and spectroscopic data, using a modified version of 
the EAZY code (Brammer, van Dokkum, & Coppi 2008). 
The slitless spectra require that the template fitting explicitly 
takes the morphologies of the galaxies into account; therefore, 
the models are convolved with the F140W direct image aver- 
aged in the spatial direction. In our fitting procedure emis- 
sion lines are "automatically" corrected for underlying stellar 
absorption. Due to the nature of slitless grism spectroscopy, 
many spectra are contaminated to some extent by overlapping 
spectra of neighbouring objects. The aXe package provides 
a quantitative estimate of the contamination as a function of 
wavelength, which can be subtracted from the spectra. 

Stellar masses were determined using the FAST code 
(Kriek et al. 2009a), using Bruzual & Chariot (2003) models 
and assuming a Chabrier (2003) stellar initial mass function 
(IMF). We find that the stellar masses differ by < 0.05 dex 
whether or not the spectra are included in the fits. Rest-frame 
U — V colors were determined directly from the best-fitting 
spectral energy distributions, following the procedures out- 
lined in, e.g., Wolf et al. (2003) and Brammer et al. (2011). 
Effective radii and Sersic (1968) indices were determined 
from the F140W images using GALFIT (Peng et al. 2002), 
using synthetic PSFs and averaging the measurements from 
each of the dithered exposures. The individual measurements 
are generally consistent to within < 0.05 dex. 
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Fig. 1. — Color-mass relation of galaxies at 1 < z < 1.5 in the 
NEWFIRM Medium Band Survey (greyscale). The 3D-HST galax- 
ies are overplotted, color-coded by their Ha equivalent width. The 
3D-HST galaxies span a similar range in color as the much larger 
NMBS sample with M > 10" M©. Star-forming galaxies are typi- 
cally bluer than galaxies with weak Ha. The line is the separation 
between blue and red galaxies from Borch et al. (2006), extrapolated 
to z = 1.25. 



3. A MASS-LIMITED SPECTROSCOPIC GALAXY 
SAMPLE AT Z= 1-1.5 

For this initial paper, we selected galaxies from the 3D- 
HST survey with stellar masses M stal - > 10 11 M Q , grism red- 
shifts 1 < z < 1.5, and low contamination (less than 20 % of 
the measured flux). The only other spectroscopic sample of 
this kind is the Gemini Deep Deep Survey, which used op- 
tical spectroscopy (Abraham et al. 2004). The upper limit 
on the redshift selection is the highest redshift for which Ha 
falls in the G141 spectral range. 14 Besides Ha, strong Mg, 
Na, and TiO absorption features fall in the observed wave- 
length range, allowing us to measure accurate redshifts for 
both star-forming and quiescent galaxies. From 68 available 
HST pointings (approximately 50 % of the full survey) this se- 
lection gives a sample of 34 galaxies. We note that our sample 
is biased against close pairs, mostly due to the strict contami- 
nation criterion. Fifteen of these objects have a spectroscopic 
redshift previously measured with a ground-based telescope. 
There is good agreement between the grism redshifts and the 
spectroscopic redshifts: the la scatter is 0.004 x (1+z). 

In Fig.Q]the rest-frame U — V colors of the galaxies in our 
sample are compared to those of 1 < z < 1.5 galaxies over 
a large mass range in the NEWFIRM Medium Band Survey 
(NMBS; Whitaker et al. 201 1). The 3D-HST galaxies span a 
similar range in rest-frame U — V color as massive galaxies in 
the NMBS. This range is much smaller at M > 10 11 M than 
at lower masses: all massive galaxies are red compared to blue 
cloud galaxies at M — 10 10 M Q , which have U-V ~ 0.8. 

14 At the resolution of the 3D-HST spectra Ha and [N II] A6583 cannot be 
distinguished. Throughout the text "Ha" refers to the combination of these 
two lines. 
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Fig. 2. — HST/WFC3 F140W images and G141 grism spectra of galaxies with M sla r > 10 11 M© and 1 < z < 1.5. The spectra were averaged 
in the spatial direction to optimize the S/N ratio, and smoothed by a boxcar filter for presentation purposes. The v-axis is in units of F\; each 
spectrum was normalized and offset with respect to the others. Squares indicate broad- and medium-band photometry from public catalogs (see 
text). The orange spectra are best-fitting EAZY models (Brammer et al. 2008) for the continuum emission and for Ha. In the grism spectral 
range the models were convolved with the morphologies of the galaxies to properly model the spectral resolution. The galaxies are ordered by 
decreasing Ha equivalent width. The spectra are of high quality. The fraction of galaxies with strong Ha emission is much higher than in the 
nearby Universe. Galaxies with strong Ha emission are often two-armed spiral galaxies. Galaxies with weak or undetected Ha typically have 
an early-type (E, SO, or Sa) morphology. 
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4. SPECTRAL FEATURES 

The WFC3/G141 spectra and WFC3/F140W images of the 
galaxies are shown in Fig. [2] ordered by decreasing rest-frame 
Ha equivalent width. We detect Ha emission in 20 galaxies, 
or 59 % of the sample. The rest-frame equivalent widths for 
the detected galaxies range from 10 A to 130 A. The immedi- 
ate implication is that a substantial fraction of massive galax- 
ies at 1 < z < 1.5 is likely forming stars. The Ha lines in 
nearby galaxies in this mass range are substantially weaker: 
the fraction of galaxies with EWn a > 10 A above the same 
mass limit at z = 0.1 is only ss 10% (Tremonti et al. 2004). 
Note that the Ha lines are typically broad in the 3D-HST 
spectra. This is not due to velocity broadening but rather to 
"morphological broadening": in slitless spectroscopy spectral 
features are effectively images of the galaxy in that particular 
wavelength. This will be discussed further in §[6] 

The galaxies with weak Ha emission have strong stellar 
absorption features typical of intermediate age stellar popu- 
lations. This is demonstrated in Fig. [3] which compares the 
averaged rest-frame spectra of galaxies with strong and weak 
Ha. The stellar absorption features Mg&, NaD, and several 
TiO bands are clearly detected in the stacked spectrum of the 
weak Ha galaxies - for the first time at these redshifts. The 
colored lines in Fig.[3]are stellar population synthesis models 
of Vazdekis et al. (2010) with metallicity [M/Z] =0.22 and 
different ages. The best-fitting age for the galaxies with low 
star formation is 2 Gyr; models with ages of up to 4 Gyr also 
provide good fits. We note that this is an average and that the 
galaxies probably have a range of ages. 

As can be seen in Fig.[T]the galaxies with strong Ha emis- 
sion are slightly bluer (by ~ 0.3 mag in U — V) than the 
galaxies with weak Ha. The fact that the color difference 
is relatively small is because the star-forming galaxies have 
more dust: gaxies with EWho < 10 A have a mean best-fit 
(Ay) = 0.5 ±0.1 whereas galaxies with EWh q > 10 A have 
(A v ) = 1.0 ±0.2 (see also, e.g., Labbe et al. 2005; Papovich 
et al. 2006; Brammer et al. 201 1). 

5. STAR FORMATION, STRUCTURE, AND THE 
DIVERSITY OF MASSIVE GALAXIES 

Besides the number of objects with strong Ha emission, a 
striking aspect of Fig.|2]is that the morphologies of the galax- 
ies correlate with the Ha line strength. Massive 1 < z < 1.5 
galaxies with the highest star formation rates tend to be "grand 
design" spiral galaxies, whereas those with no detected Ha 
emission tend to be early-type galaxies. This result is quali- 
tatively similar to trends at z = 0, although massive galaxies 
with EWn a > 10 A are rare in the nearby Universe (see § 3). 

We quantify this correlation between galaxy structure and 
star formation rate in the top panel of Fig. U which shows 
the relation between EWh q and Sersic (1968) index. The Ser- 
sic index is a quantitative measure of galaxy structure, and 
is a proxy for the bulge-to-disk ratio: galaxies dominated by 
disks have n ~ 1 and galaxies dominated by bulges have n ~ 4. 
Splitting the galaxies in two equal bins, the median Sersic in- 
dex of the 17 galaxies with EWh q > 12 A is 2.2, whereas it is 
4.2 for the 17 galaxies with EWh q < 12 A. According to the 
Mann- Whitney test this difference is significant at the > 99 % 
confidence level. We infer that star formation in massive 
galaxies at 1 < z < 1 .5 typically takes place in disk-dominated 
(spiral) galaxies. This result is consistent with recent sugges- 
tions that many star-forming galaxies at z > 1 are "scaled-up" 
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Fig. 3. — Stacked spectra of the 20 massive galaxies at 1 < 
z < 1.5 with EWhq > 10 A (top panel) and of the 14 galaxies with 
EWhq < 10 A (bottom panel). The spectra were normalized by di- 
viding them by a second-order polynomial fit. Colored lines are 
stellar population synthesis models with different ages (see text), 
smoothed to the resolution of the stacked spectra. The strong stellar 
absorption features of the galaxies with weak Ha require luminosity- 
weighted ages of 2 -4 Gyr. 

versions of nearby galaxies (e.g., Elbaz et al. 201 1). 

The bottom panel of Fig. U shows the relation between 
EWhq and the inferred velocity dispersions of the galaxies, 
which are derived from their stellar masses, effective radii, 
and Sersic indices following Bezanson et al. (2011). The in- 
ferred dispersions are a measure of the compactness of the 
stellar components of galaxies as they are proportional to 
\jMjr e (with a Sersic-dependent correction factor). Again, 
there is a clear relation between these quantities: the median 
inferred dispersion of the half of the sample with the strongest 
Ha is 174^25 km/s whereas it is 299^3 km/s for the galaxies 
with the weakest Ha. The relation between EWh q and dis- 
persion has much smaller scatter than the relations of EWu a 
with effective radius and mass separately. The trend in Fig. [4] 
is qualitatively consistent with the relation between estimated 
specific star formation rate (SSFR) and inferred dispersion 
found by Franx et al. (2008). 

Figure 2] showcases the large variety of galaxies. The range 
in EWhq among massive galaxies with detected Ha is a fac- 
tor of 12, and it is obviously even larger when undetected 
galaxies are taken into account (see Fig. [3). The structure 
and morphologies of galaxies also show a large range, going 
from large spiral galaxies with high EWh q to compact early- 
type galaxies with low EWh q . The Sersic indices span the 
full range between disk- to bulge-dominated values (« ~ 1 to 
n ~ 4) and the inferred dispersions range from < 150 km s" 1 to 
> 400^8-'. 
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Fig. 4. — Top panel: Relation between EWh , Sersic index, and morphology. The F140W images are color-coded according to their average 
U-V color. Bottom panel: Relation between EWho,, velocity dispersion, and morphology. The velocity dispersions were inferred from 
the stellar masses, effective radii, and Sersic n values of the galaxies. The errorbars denote the typical uncertainties, as determined from a 
combination of formal errors and varying the fitting methodology. Sersic index, velocity dispersion, and morphology all correlate with the 
strength of Ha. All these parameters span large ranges: massive galaxies exhibit a striking diversity at 1 < z < 1.5, ranging from large 
star-forming spiral galaxies to compact quiescent galaxies. 
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6. DISCUSSION 

We have shown that the 3D-HST Treasury program is pro- 
viding high quality near-IR spectroscopy and accompanying 
imaging for galaxies in the crucial epoch 1 < z < 3, when the 
cosmic star formation rate peaked (e.g., Bouwens et al. 2007). 
In this first paper we focused on a stellar mass-limited sample 
of 34 galaxies at z = 1 - 1 .5, to study the properties of massive 
galaxies when the Universe was 4-6 Gyr old. This is the first 
study of the Ha emission and rest-frame optical morphologies 
of a complete, mass-limited galaxy sample at z > 1 . 

The most striking result of our study is the diversity of the 
spectra and structure of massive galaxies at z > 1: a large 
fraction has strong Ha emission, whereas others have ab- 
sorption features characteristic of relatively old stellar pop- 
ulations. Similarly, the morphologies, Sersic indices, and im- 
plied velocity dispersions show a large range. These results 
are broadly consistent with previous studies that were based 
on photometric redshifts, the SED shapes of galaxies, and/or 
imaging of lower quality (e.g., Franx et al. 2008; Williams 
etal. 2010; Wuyts et al. 2011b; Weinzirl et al. 2011). The 
simplest interpretation is that at z > 1 we are entering the 
epoch when massive galaxies were undergoing rapid evolu- 
tion. Specifically, the star-forming disks may be progenitors 
of some of today's most massive SO and Sa galaxies. 

Noeske et al. (2007) suggested that the star formation rates 
of galaxies are tightly coupled to their mass and redshift, and 
it is interesting to compare the range in EWh q in our sample 
to the scatter in their "star formation main sequence". At fixed 
mass Noeske et al. (2007) find a 68 % range in star formation 
rates of a factor of ~ 4 (among galaxies with clear signs of star 
formation); if we limit the sample to the central 68 % of the 
distribution of galaxies with EWh q > 10 A we find a larger 
range of a factor of six. More to the point, within our sample 
the scatter in EWh q can be significantly reduced by consid- 
ering the structure of the galaxies. The 68 % range in EWh q 
among the 11 galaxies with inferred a < 200 km s" 1 is only a 
factor of 1.7. We therefore follow earlier work in suggesting 
that velocity dispersion (or surface density) is a more funda- 
mental parameter than mass in determining the properties of 
galaxies (see, e.g., Kauffmann et al. 2003b; Franx et al. 2008; 
Bezanson et al. 201 1). 

The growth rate of the star-forming galaxies is substantial; 
using standard prescriptions to correct for extinction toward 
Hll regions (Calzetti et al. 2000; Wuyts et al. 2011a), we 
find a median stellar mass increase due to star formation of 
~ 50 % per Gyr for the galaxies with EWh q > 10 A. An im- 
portant question is where in the galaxies the star formation 
is occurring, that is, which structural component of massive 
galaxies is in the process of formation at 1 < z < 1.5. Due 
to the nature of grism spectroscopy the 3D-HST data provide 
2D emission line maps at the spatial resolution of HST Two 
examples are shown in Fig. [5] in these galaxies the star forma- 
tion appears to trace the spiral arms, similar to spiral galaxies 
in the nearby Universe. A quantitative analysis of the spatial 
extent of the emission line gas is beyond the scope of this Let- 
ter, but we note here that in cases such as those shown in Fig. 
[5] the spatial extent of the Ha emission rules out dominant 
contributions from active nuclei to the integrated line fluxes. 

A key open question is what drives the diversity of mas- 
sive galaxies at z > 1 . At fixed stellar mass, we see large, star 
forming spiral galaxies and very compact galaxies in which 
star formation has apparently ceased. If AGN feedback is re- 



sponsible for shutting off star formation in massive galaxies 
it is clearly more effective in some galaxies than in others. It 
may be that AGN feedback correlates with black hole mass, 
which correlates better with velocity dispersion than with stel- 
lar mass (e.g., Magorrian et al. 1998). It will also be interest- 
ing to study correlations with other parameters, such as the 
environment, at fixed stellar mass and at fixed (inferred) ve- 
locity dispersion. Finally, it will be important to extend this 
study to lower masses and to higher redshifts. Star forming 
galaxies that have been studied at z > 2 tend to have higher 
EWhq and also more irregular morphologies than the galax- 
ies studied here (e.g., Erb et al. 2006; Kriek et al. 2009b; 
Forster Schreiber et al. 201 1), and it will be interesting to see 
whether a broader selection of galaxies would only increase 
the dynamic range in Fig.|4]or (also) increase the scatter. 

Such studies will be possible with the full 3D-HST dataset, 
which will provide WFC3 and ACS spectroscopy and imag- 
ing of large samples of galaxies at z > 1. Combined with 
deep imaging from the CANDELS project and the wide array 
of ancillary data in the survey fields, 3D-HST will provide a 
qualitatively new way of surveying the heyday of galaxy for- 
mation at 1 < z < 3. 
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Fig. 5. — Spatial distribution of continuum emission in the F140W 
filter compared to the spatial distribution of the Ha line emitting gas. 
The Ha images were made by subtracting a polynomial fit from the 
two-dimensional G141 grism spectra. The boxes are 4."7 x 4"7, or 
39 kpc x 39 kpc. The line profiles show the flux distributions aver- 
aged along the y-axis. The line emission in these particular galaxies 
is extended, and the galaxies are building up their disks. We also find 
galaxies with more compact Ha distributions. 
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